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We discuss two new observational techniques that use observations of the 
J2 I Cosmic Microwave Background (CMB) to place constraints upon the mass, dis- 

tance, and size distribution of small objects in the Kuiper Belt and inner Oort 
I Cloud, collectively known as Trans- Neptunian Objects (TNOs). The first new 

technique considers the spectral distortion of the isotropic, or monopole, CMB 
by TNOs that have been heated by solar radiation to temperatures above that 
Q I of the CMB. We apply this technique to the spectral measurements of the CMB 

^ I by the Far Infrared Absolute Spectrophotometer (FIRAS) on the Cosmic Back- 

ed I ground Explorer (COBE). The second technique utilizes the change in amplitude 

of the TNO signal due to the orbital motion of the observer to separate the TNO 
^ , signal from the invariant extra-galactic CMB and construct a map of the mass 

' distribution in the outer Solar System. We estimate the ability of future CMB 
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, , . 1. Introduction 

■ While we know that in some manner dynamical instabilities in the proto-planetary disk 

led to the format ion of planetesimals and planets, the precise detai ls of this process are not 



experiments to create such a map. 



well understood (iGoldreich et al.l 12004 iLissauer fc Stevensonll2005l ). Luckily, just as nature 



has provided us with the Cosmic Microwave Background (CMB) as a relic of the early 
universe, there are similar relics of this chaotic period of formation. Objects in the outer 
Solar System, collectively known as Trans-Neptunian Objects (TNOs), with semi-major axes 



^Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138; E-mails: 
dbabich , cblake , est einha@cf a. harvard . edu 

^Cahfornia Institute of Technology, Theoretical Astrophysics, MC 130-33, Pasadena, CA 91125 

•^Harvard Origins of Life Initiative Fellow 



- 2 - 



40 AU < a < 10^ AU are believed to provide clues that may help us understand the early 
history of the Solar System. The number density and mass distribution of TNOs contain 
valuable information about the properties of the proto-planetary disk out of which the Solar 
System formed and the dynamical properties of these TNOs can constrain migratory motion 
of the planets. We show in this paper that the CMB itself may be the key to unlocking the 
mysteries of these objects. 

To date, the entirety of our observational evidence about the structure of the Solar 
System beyond Neptune consists of observations of comets that have been perturbed into 
the inner Solar System and of r elatively large, nearby objects detected through their reflected 
sunlight (IBernstein et al.ll2004j ). Because of the strong dependence of the reflected sunlight on 
the object's distance and radius (oc cx R~^ ), it is very difficul t to find objects smaller or 



more distant than, for example, (90377) Sed na (IBrown et al 



200^, a body 1180 -1800 km in 



diameter detected at 90 AU, or 2003 UB313 jBrown et al.lboosh . a body 2400 km in diameter 



detected near its aphelion distance of 96 AU. 

Long Period Comets (LPC) associated with the Oort Cloud originate in its outermost 
reaches because objects with large semi-major axes can be preferentially ejected into the 
inner Solar Systeno. The observed aphelia of L PCs place t he outer limits of the Oor t 



Cloud at a distance of approximately 25000 AU (lOort 



I950I : iMarsden fc Sekaninal Il971f ) 



while simulations of the Oort Cloud's formation, assuming that the Sun formed in a star 
cluster, suggest that a significa r it amount of ma s s may lie in an in ner Oort Cloud, which is 
located at 1000 AU JHillslll98ll : lFernandez|[l997l : IPones et al1l2004h . To this point, we have 
no methods capable of detecting smaller objects in the Kuiper Belt or any objects in the 
inner Oort Cloud. 

In this paper, we develop two new methods capable of exploring the distribution of 
TNOs. Both utilize the CMB as a standard spectral template, namely a blackbody with a 
given temperature emanating from a surface at infinity. The Sun heats Kuiper Belt and Oort 
Cloud objects to a temperature above that of the CMB. These objects extinguish part of the 
CMB and emit blackbody radiation, inducing a spectral distortion in the observed CMB. 
By constraining these spectral distortions we can constrain the radial and mass distribution 
of trans-Neptunian objects. 

Some similar topics have previously been discussed. Detections of zodiacal light, both 
through its refiected optical sunlight and its emitted infrared radiation, have led to detailed 



^Perturbations by Jupiter are large enough to cither eject the new LPC on a hyperbohc orbit or keep it 
within 1000 AU of the Sun. This originaUy led Oort to postulate that comets on orbits with large aphelia 
are always entering the inner Solar System for the first time. 
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models of inter-planetary dust in the Solar System (IFixsen fc Dwekl |2002| ). This dust is 
believed to come from the comae of comets and the collisional debris of asteroids. The 
infrared emission of Kuiper Belt objects h as been employed in a sinailar fashion to this 



paper in order to con strain the Kuiper Belt (IKenyon &: Windhorstll200ll : iTeplitz et al.lll999 



Backman et al.lll995l ).The main advantage of using the CMB as the source of background 
radiation is that its properties are much better understood and more uniform than the 
far-infrared background, allowing for a greater precision. Additionally, the longer relevant 
wavelengths of the CMB spectrum enable us to apply the same test to the inner Oort Cloud, 
which is at too low of a temperature to emit any substantial amount of radiation in the 
infrared bands. 

In §[2|, we calculate the blackbody radiation emitted by objects in the outer Solar System 
including a discussion of the temperature of distant objects. In § [31 we use FIRAS to place 
constraints on the total mass and size distribution of the smallest objects in the Kuiper belt. 
In § m we propose a new method of using WMAP to develop a two-dimensional map of the 
smallest objects in the Kuiper Belt, holding out the possibility that for sufficient mass we 
may be able to determine the distance in addition to the mass and direction. 



2. Outer Solar System Spectral Distortions 



The extremely high interaction rates between matter and radiation during the early 
universe cause the primordial CMB spectrum to have a blackbody distribution to incred- 
ibly high precision (IMather et al.l Il992l : iPeebled 1 19931 ). However, as the universe expands 
and these interaction rates decrease, new sources of energy injection can distort the CMB 
spectrum and photon-baryon plasma is now incapable of producing a blackbody spectrum 
again. In addition to energy injection at high redshifts, there are several low r e dshift pro- 
cesses that are capabl e of distorting the CMB spectrum (see iTegmark fc Silk! (jl99J) and 
Fixsen fc Kogutl (|2002| ) for an overview). In this paper we will try to constrain the proper- 
ties of TNOs from the absence of spectral distortions, so any additional process that may 
produce these distortions does not affect our ability to place these constraints on the outer 
solar System, but in principle strengthens our conclusions. We will consequently ignore these 
other effects in this paper. 

It should be emphasized that these spectral distortions are not the sam e as the chemical 
potential (fj,) or Compto n-?/ distortions often discussed in the literature (ITegmark fc Silk 



1994 : iFixsen et al.l 119971 ). Of course, there is a level of degeneracy between the spectral 



distortion considered in this paper, namely a weighted sum of blackbodies, and the chemical 
potential and Compton-y distortions, so our signal will result in a non-zero chemical potential 
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and Compton-?/ distortion if either of these models was assumed in the analysis. 



2.1. Calculation 

We are interested in the aggregate blackbody emission from objects at a given distance 
in the outer Solar System. The distorted CMB intensitjQ spectrum due to these objects will 
be 

5I, = t[B,{Ttno)- B,{Tcmb)], (1) 

where the optical depth, or equivalently the geometric covering fraction since we are working 
in the geometric optics limit, is expressed as 

r = y ^n(M)rfM, (2) 

where R is the radius of the object and D is the approximately identical distance to the 
TNO from both the Sun and the observer. These two distances are not exactly equal and 
the TNO-observer distance may be time-variable. We will explore this possibility in ^ 

In this paper we assume that all of the TNOs are at a single, unknown distance which is 
then constrained by observations. Of course, TNOs can have quite different eccentricities and 
semi-major axes so at a given time there will be a radial distribution of TNOs. Since both the 
geometric covering fraction and TNO temperature decrease with increasing distance from 
the Sun, we will be most sensitive to the closest TNOs. This helps to justify our assumption 
of a single distance. Moreover, the TNO temperature varies with distance from the Sun 
so the radial distribution of TNOs will produce a more complicate spectral distortion than 
implied by Eq. ([1]). With good enough data this radial distribution could be constrained by 
constraining the distribution of TNO temperatures. However, with the limited sensitivity 
of current data our assumption of a single distance is sufficiently accurate. In §1] we will 
slightly relax this assumption by allowing the TNO distances to be anisotropic; therefore, 
all of the TNOs in a given pixel will be at the same distance, which will be different than 
the TNO distance in other pixels. 



^ There should be no measurable polarization signal as the mean CMB and the TNO emission are both 
blackbody radiation. Any polarization signal produced by features in the surface of the TNOs will vanish 
because the experiment's beam will contain many TNOs that are randomly oriented with respect to one 
another. If the TNOs are aspherical and have a large mag netic suscepti bility it is possible that they may 
align with the magnetic field via the Greenstein-Davis effect (jPraind 120031 ). An analysis of this possibility is 
beyond the scope of this paper. 
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The size distribution of TNOs is specified by their mass function, n(M), with total mass 

Mtotai = / n{M)MdM. (3) 



For Kuiper Belt objects, the heating is almost entirely solar, so we can calculate the 
equilibrium temperature of these objects. In the case of Oort Cloud objects, their equilibrium 
temperature even at large distances is higher than that of the CMB, even though at 25000 
AU the solar contribution is tiny. Interstellar processes, as well as the absorption of CMB 
photons, results in a floor for the TNO temperature. The equilib rium temperature of Oort 



Cloud objec ts is assumed not to drop below 5-6 K, as claimed by iMumma et all (119931 ) and 



SternI (120031 ) . In this work we will calculate the temperature of TNOs by assuming thermal 
equilibrium with the Sun. As mentioned above, this assumption breaks down for TNOs with 
extremely large semi-major axes (a > 10000 AU), but objects at such a large distance will 
produce a signal well below current instrumental sensitivities and therefore are irrelevant for 
our work. Assuming thermal equilibrium with the Sun, the temperature of a TNO is 



TNO 



R 



4^2 



1/4 



^0, (4) 



where Rq is the Sun's radius, Tq the temperature of the Sun's photosphere, and A the TNO's 
albedo. At distances appropriate for the Kuiper Belt, D = 40 AU, we find a temperature of 
T = 43 K and at inner Oort Cloud distances, D = 1000 AU, the temperature is T = 9 K. 

The total distortion will depend on the following parameters of the Oort Cloud and 
Kuiper Belt. Here we describe our choices for these parameters, as well as theoretical and 
observational constraints. 

• Density - We assume a density p = 1 g cm~^. This number depends on the porosity of 
the TNOs and may vary by a factor of a few. 

• Albedo - We will assume an albedo of A = 4% as has been measured for comets and is 



relevant for dirty ice (ILuu fc Jewittll2002l ). The constraints only depend on the albedo 



as (1 — Ay/^ so this uncertainty does not strongly affect our results. 
Mass Function - We will assume a broken power 



■^Typically the differential mass function is expressed in terms of comet radius, not mass. The power law 
with respect to radius is related to the values used in this paper as amass = (^radius + 2)/3. The canonical 
faint-end power-law exponent of aradius — 3.5 implies amass — 1.83. 
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n{M) 
n{M) 



AM, 



AM-'', Mmin < M < Mbr 

-f+^M-f, M^,<M<M^,^. 



br 



The appropriate power law exponents are not well known. We will find that our results 
are quite sensitive to a, the low mass slope, because for a given total amount of mass in 
either the Oort Cloud or Kuiper Belt a steeper mass function will increase the geometric 
covering fraction. It is difficult to observationally constrain the mass function of comets 
due to the uncertainties in the dynamics of comet's comae. Collisional equilibrium, 
along with the assumption that the strength of the object is independent of size, would 
lead us to expect a = 11/6 (jPan fc Sarill2005l ). If the larger objects are more difficult to 
shatter, as would occur if the gravitationally binding energy were very important, then 
a < 11/6. Since there is so much uncertainty in a we will present results for several 
values. The value of /3 has little effect upon our final results if /3 > 2. We will adopt 
a single value of = 1 3/6, which is con sistent with the current observational results 
( iBernstein et al.l |2004J : iPan fc Saril 120051 ) . In ^ we will discuss how our conclusions 
depend on this assumption. 

Lower Mass - The minimum size of surviving objects is determined by the ability of 
Robertson-Poynting drag to eliminate the smallest objects as we ll as the propertie s 



of collisions that can fragment larger objects into smaller ones (iBurns et al. 



19791). 



We will take this lower mass to correspond to objects of radius 1 mm. This quantity 
depends on the detailed formation history of the Kuiper Belt and the Oort Cloud and 
is difficult to calculate. This quantity is partially degenerate with the total mass of 
the outer Solar System, which is also unknown, so for the purposes of this paper we 
will fix it at 1 mm. 

However, we should note that the efficiency of blackbody emission is suppressed when 
the radiation wavel ength is larger th an the size of the object emitting the radiation due 
to Kirchoff's Law (iGreenbergl Il978l ) . This suppression of the TNO emissivity, which 
will change the low frequency portion of the CMB spectrum, may be used to constrain 
the minimum mass. We will discuss this in more detail in the ^ where we mention a 
new proposed CMB experiment to measure these low frequencies. 

Upper Mass - We assume M^ax = 10^^ x Msarth = 6 x 10^^ g as the upper mass. Our 
results are insensitive to changes of even many orders of magnitude in the upper mass 
limit for (3 > 2. If /3 < 2, a. significant fraction of the mass would be at the high-mass 
end and our limits on the low-mass end would not directly translate to limits on the 
overall mass in the outer solar system. 
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Break Mass - Mbr = 3.2 x 10^^ g (corresponding to i?br = 2 x 10^ cm). The break mass 
is analytically calculated by determining the large st TNOs that can be in collisional 
equilibrium over the lifetime of the Solar System (jPan &: Sarill2005l ). Like the power 
law indices, this parameter is also uncertain but, like the upper mass, our results are 
relatively insensitive to its value. Observations with the Hubble Space Telescope have 
constrained the break radius to be -Rbr < 20km (IBernstein et al.l 120041 ). 



Total Mass - The total mass Mtot is one of the things which we seek to constrain with 
this calculation. Theoreti cally the total mass should be between (0.1 — 100) x Msarth (cf. 
Stern fc WeissmanI (l200ll )). depending on the surface mass density of proto-planetary 
disk and the details of the growth of the planets. 

Distance - We consider what limits can be set for distances from D = 40 AU to D = 
10000 AU. The strength of the effect roughly scales as so we are only sensitive to 
the inner Oort Cloud and not the outer Oort Cloud from where LPC are observed to 
originate. 



2.2. Strategies 

In this subsection we will outline the two strategies that can be used to constrain the 
mass distribution in the outer Solar System via their spectral distortions to the CMB. For 
small optical depth r, the observed intensity fluctuation in the CMB towards a given direction 
h is 



= [1 - T{h)] [B,{f CMB + AT(n))] (5) 

+n-V + T{n)By{Toort{n)) + iVinstr, 



corresponding to the extragalactic CMB, including the temperature anisotropy (AT(ri))), 
extinguished by the TNOs in the beam; the Doppler effect, assuming that the observer has 
velocity v with respect to the CMB restframe; TNO blackbody emission; and instrument 
noise, which is statistically stationary but random. In what follows, we will assume that 
the anisotropy produced via the Doppler effect, while in principle large, is calculable and 
therefore does not affect our ability to extract or constrain the TNO component of the signal. 
Subtracting off the mean CMB blackbody, the intensity fluctuation in a given direction is 



SL(n 



V 

X 



= AT(r2) + r(n) (6) 

{Bu(Toort{n)) - B,{TcMB + AT(n))] + N^^.^,, 
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Unfortunately for our purposes, the CMB has small statistical temperature fluctuation^l on 
the order of 1 part in 10^. We would like to constrain the TNO optical depth, T(n), in a given 
direction. This is complicated by the random terms in Eq. corresponding to instrument 
noise and CMB temperature anisotropies. 

We should be able to directly distinguish the > 5 K blackbody spectrum from the CMB 
blackbody spectrum, especially when we include data at higher frequencies {u > 200 GHz), 
at a level set by either the CMB temperature anisotropies or instrument noise. This would 
enable us to directly make a map of the TNO mass density. If the instrument pixel noise 
dominates both the intrinsic variation due to the CMB temperature anisotropies and signal 
from the outer Solar System, then the constructed TNO map would have a signal-to-noise 
ratio of less than one and would be useless. In this case, we can reduce the instrument noise 
by smoothing the CMB map and averaging together neighboring pixels; this will increase the 
pixel signal-to-noise ratio. The extreme limit of the procedure is to calculate the mean CMB 
spectrum. If the CMB temperature anisotropies dominate the instrument noise it is possible 
to reduce, but not completely eliminate, the variance of the temperature anisotropies. This 
technique is based on differencing the observations of a given position of the Celestial sphere 
made at different locations within the Solar System. It will be discussed in detail in §H 
Using this technique it may be possible to make a two- or even three-, dimensional map of 
the distribution of TNOs. 



3. Isotropic Spectral Distortion 

3.1. Observational Constraints 

We use data from the Far Infr ared Absolute Spect rophotometer (FIRAS) on the Cos- 
mic Background Explorer (COBE) (IFixsen et al.l 119971 ) to place limits on deviations from 
a blackbody spectrum due to TNOs. FIRAS measured the mean, or monopole, CMB fre- 
quency spectrum so we will average the signal originating in the outer Solar System over 
the entire sky in order to compare with their results. The Oort Cloud is assumed to be 
isotropic since the observed distribu tion of inclination angles of n ew LPCs is approximately 
uniform (IMarsden fc Sekaninalll97ll), although the Galactic tide (IHeisler fc Tremaindll986l ) 
and stellar perturbations (IBabich fc Loebl 120071 ) should make it aspherical. 



"^This implies that it is impossible to predict the exact pattern of observed temperature anisotropies. We 
can only predict the two-point correlation function (the scale-dependent variance) which parameterize the 
probability distribution function of which the observed CMB is a random realization. 
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The Kuiper Belt is highly anisotropic, since it is primarily located near the ecliptic 
plane. The FIRAS results we use have been averaged over the full-sky in each frequency 
band and therefore the signal from the Kuiper Belt must also be averaged over the full sky. 
Since we know that the Kuiper Belt is located at small inclination, it is less optimal to use 
the full sky averaged FIRAS spectrum. This will still allow us to place constraints on the 
Kuiper Belt, even if they are not the best possible. 

Since FIRAS did not detect any spectral distortion in the CMB we can only place upper 
limits on the basic parameters of the outer Solar System. There are a multitude of processes 
that can produce spectral distortions but we will ignore them as we are simply using the lack 
of detection to constrain the properties of the outer Solar System. Note that the inclusion 
of these other effects would only make our upper limits stronger. 

Figures [T] & [2] display constraints on the total mas^ and the distance of the Kuiper Belt 
and inner Oort Cloud, respectively, from a non-detection of CMB spectral distortions in the 
FIRAS data. The confidence contours are determined by calculating the change in the 
of the FIRAS data when a component originating from the o uter Solar System is included. 



The errors used to calculate the are taken from Table 4 of (IFixsen et al.lll996[ ) for the low 
frequency FIRAS data which extends up to 630 GHz. 

The curves are shown for several values of a since the constraints strongly depend on 
the low mass end of the planetesimal mass function. Even though the Kuiper Belt primarily 
lies in the ecliptic plane and the FIRAS data is the mean CMB spectrum, averaged over the 
entire sky, we can still use this data to constrain, in a sub-optimal manner, the properties 
of the Kuiper Belt. These constraints correspond to 95% confidence limits. Also shown is 
a dynamical const raint on the tota l mass in the Kuip er Belt from constraints on the orbit 



of Halley's comet (IHogg et al.lll99ll : iHamid et al.lll968l ). These figures demonstrate that our 



technique is competitive with dynamical limits on the Kuiper Belt mass and may place more 
stringent constraints if the slope of the low mass end of the mass function is steep enough. 
Hence, if amass = 3.5 {aradius = 1-83) then the dynamical CMB results place limits to the 



^Our technique is mainly sensitive to low mass objects so we must be very cautious when interpreting 
these results. First of all, we are assuming that the mass function is well modeled as a broken power-law 
over a very wide range in mass. Second, Poisson fluctuations become important at the high mass end, where 
the expected number of objects, as extrapolated from the mass function, approaches unity. Additionally, 
run-away growth during the core accretion phase of formation can grow a single object to a much larger mass 
than any other object in its vicinity. A simple extrapolation of the mass function will miss these objects. 
These last two effects can cause the actual mass of the Kuiper Belt or Oort Cloud to differ significantly from 
our constraints. Fortunately, these high mass objects are more easily detected in optical surveys and so our 
technique is complimentary to existing methods. 
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TNO mass that are as tight as the dynamical hmits from Halley's comet's orbit, and exclude 
M > IMq in general. This technique may be the only way we can constrain the properties 
of the inner Oort Cloud. Of course, our technique is only sensitive to the low mass end of 
the distribution and so care must be taken when interpreting these results. 



2 Q-4 I I I I I I I I I I I I I 

40 60 80 100 

Distance (in AU) 

Fig. 1. — Excluded regions of total mass and distance for the Kuiper Belt for different mass 
function slopes a. The curves correspond to a = 2 (black, solid); a = 1.9 (blue, dashed); 
a = 1.8 (red, long-dashed); a = 1.667 (black, dot-dashed); a = 1.5 (blue, dotted) and 
a = 1.333 (red, dashed). The region above the curves is excluded at 95% by the COBE 
FIRAS data. Also shown is a dynamical limit determined by constraints on perturbations to 
the orbit of Halley's comet (heavy black, solid curve) and, for reference, a line corresponding 
to a Jupiter mass (heavy red, horizontal line). The power law with respect to radius is 
related to the values used in this paper as amass = {(^radius + 2)/3. 
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Fig. 2. — Excluded regions of total mass and distance for the inner Oort Cloud for different 
mass function slopes a. The curves correspond to a = 2 (black, solid); a = 1.9 (blue, 
dashed); a = 1.8 (red, long-dashed); a = 1.667 (black, dot-dashed); a = 1.5 (blue, dotted) 
and a = 1.333 (red, dashed). The region to the upper left of the curves is excluded at 95% 
by the COBE FIRAS data. For reference we also include lines corresponding to a solar mass 
(heavy black, horizontal line) and a Jupiter mass (heavy red, horizontal line). 
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3.2. Low Frequency Observations 



A new experiment, tlie Absolute Radiometer for Cosmology Astrophysics and Diffuse 
Emiss ion (ARCADE), li as been proposed to constrain low frequency CMB spectral distor- 
tions (iKogut et al.ll2006l ). This experiment would measure the CMB spectrum over the full 
sky in narrow bands between 3 GHz and 30 GHz. The addition of these lower frequencies 
to the COBE FIRAS data would significantly help in constraining the mean Compton-y 
parameter and chemical potential. If the TNOs emit as perfect blackbodies, the low fre- 
quency data would not help us constrain objects in the outer Solar System. At these low 
frequencies, where both the TNO emission spectrum and the CMB spectrum are both in the 
Rayleigh- Jeans limit we would simply measure a spectrum also in the Rayleigh- Jeans limit 
at the composite temperature 



T 



CMB 



TNO- 



(7) 



Since the CMB temperature is a priori unknown it is not possible to exclude the contamina- 
tion of another blackbody at a different temperature using only low frequency observations. 
Measurements in the Wien portion of the spectrum allow us to constrain the properties of 
the outer Solar System since at these relatively high frequencies the CMB and TNO spectra 
are significantly different. However, for small enough par ticles, or low eno ugh frequencies, 
the TNO emissivity will deviate from a prefect blackbody (jGreenberglll978l ) and the emitted 
spectrum will not be well described by the Rayleigh- Jeans formula and the low frequency 
data will become useful. 

The emissivit y, of a body in therrn al equilibrium at temperature T is determined 
by Kirchoff's law (IRybicki fc Lightmanlll979l ) 



j{u)=a{v)B,{T), (8) 

where oi{v) is the absorptivity and Bi,(T) is the Planck function. In the geometric optics 
limit, which we employed in the previous subsection, the absorptivity is simply determined 
by the frequency- independent geometric cross section. When the particle size becomes com- 
parable to the radiation's wavelength and the geometri c optics limit becomes invalid the 
absorptivity becomes frequency dependent (ISpitzerlll978l ). These deviations may be used to 
constrain the minimum mass of TNOs. 



4. Anisotropic Distance Modulation 



In the previous section we considered the spectral distortion of the mean CMB when 
averaged over a long time period. Implicit in that analysis was the assumption that the 



-14- 



Sun-TNO distance (D) and observer-TNO distance [D) are identical and constant. The 
Sun-TNO distance, D, which determines the TNO's temperature, is constant during the 
hfetime of observation since the relevant periods are approximately T ~ 300 yrs for the 
Kuiper Belt and T ~ 31000 yrs for the inner Oort Cloud. D would change during the course 
of the observations if the satellite looks at the same point on the Celestial sphere from 
different locations within the Solar System. This depends on the detailed scan pattern of 
the relevant experiment. WMAP does view the same position on the Celestial sphere from 
different positions but Planck will not. 

If the scan pattern of the relevant experiment is such that it does observe the same 
point of the Celestial sphere from multiple distances, then the time-independent extragalactic 
CMB signal can be partially removed and the component arising from the TNO can be better 
constrained. In this section we study this technique by considering the various time varying 
signals and calculating potential constraints that could be produced. 



4.1. Distance-Dependent Signal 

The observed intensity fluctuation in the CMB towards a given direction h, Eq.Q, is 
dependent on the TNO-observer distance as 

/,(n, D) = [l- r{n, D)] [5,(fcMB + AT(n))] (9) 
+n ■ v{D) + T{n, D)B,{TTNo{n), D) + Ni^,,,. 

Attempting to pick out the Oort Cloud blackbody directly is difficult due to the statis- 
tical nature of the temperature anisotropics, especially if we are restricted to low frequency 
data. However, by looking at the same position on the Celestial sphere from different points 
in the CMB experiment's orbit, we can alter the distance between the observer and the TNO. 
This difference is 

- h{n,D,) = [T{n,D-:)-r{h,Dj)] 

X [B,{TTNo{n),D)-B,{f CMB + ^T{n))] 

+ h.[v{D,)-v{D,)]+Nl,,,-NU. (10) 

By taking this difference we have reduced the importance of the portion of the signal (AT x 
dB^(TcMB) / dT) that is statistically unknown. Other terms are assumed to be known to 
such a precision that they be cleanly removed. So, we consider each change in the intensity 
fluctuation accompanying a change in the CMB, as our detection limits for the TNO are set 
by the largest accompanying change. 
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4.1.1. Constancy of the CMB 

The CMB anisotropies as measured from two different locations are not necessarily 
identical. There are three effects that potentially could cause small changes: (i) each set 
of anisotropies have slightly different surfaces of last scattering and therefore depend on 
slightly different initial curvature perturbations; (ii) the path difference between the two 
different observation points can induce phase shifts in the anisotropy Fourier coefficients, 
aim, due to free streaming of the anisotropies and (in) different ordinary Sachs- Wolfe effect 
due to different values of the solar system gravitational potential at these different points. 
The first two effects produce very small changes and are negligible. Note that the ordinary 
Sachs- Wolfe effect produced by the Sun can be large: 

AT 2GMn, o 

_ w _ 1 on 1 n-8 



fcMB C2(1AU) 



1.89 X 10"^ (11) 



however, the uncertainty introduced by the effect is suppressed by the fractional uncertainty 
in the distance between the CMB experiment and the Sun. The constancy of the CMB 
anisotropies only allows us to reduce the influence of the anisotropies and their statistical 
uncertainty from 0{/ST) to (9(rAT), not completely eliminate them. 



4.1.2. Doppler Effect 



The relative velocity of the observer with respect to the CMB rest frame induces 
anisotropies into the observed CMB. For small velocities the Doppler effect can be decom- 
posed into contributions from the motion of the Solar System with respect to the CMB rest 
frame and the motion of the observer around the Suii^. The first velocity, the motion of the 
Sun with respect to the CMB restframe, will be constant on the time scale of observations. 
As we move to a different position in the Earth's orbit the Doppler contribution n ■ v will 
also change and the two Doppler effect terms in Eq. (JTOl) will not in general cancel. Even 
though the anisotropy produced by the Doppler effect is large, its uncertainty is not. The 
uncertainty in the Earth's orbit is of order 10~^^ in position a nd velocity and the uncertainty 
in the Doppler effect contribution will be of the same order (jStandishI 120041 ). 



^The motion of the observer around the Sun wiU also induce a Doppler shift in the radiation emitted by 
the TNOs. This effect is smaller than the Doppler shift of the extragalactic CMB by a factor of the TNO 
optical depth, t. 
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4.1.3. Trans- Neptunian Objects 

The temperature of a given TNO is unknown (since the distance to the object is un- 
known), but is assumed to be constant during the period of CMB observations. The fraction 
of our beam filled by TNOs depends on the variable observer-TNO distance. For all reason- 
able choices of Kuiper Belt and inner Oort Cloud parameters the beam is sparsely filled and 
therefore the TNO filling fraction will vary as 

Some fraction of the TNOs in our beam might be shifted out of the beam due to the 
parallax effect (especially when we are considering the much closer Kuiper Belt) as the CMB 
experiment changes its position within the Solar System. The fractional change of the beam 
filling fraction due to the parallax is 

where 6b is the beam width and A6 = 2AU/D is the maximum parallax angle in the 
small angle approximation. The need to minimize this fractional change will determine the 
resolution of the TNO map that we can ultimately reconstruct. By requiring this fractional 
change to be less than 1% we find 

1 ATI 

Ob > 200-^. (13) 

For the Inner Oort Cloud {D ^ 10'^ AU) 6b > 10' and for the Kuiper Belt much larger. 

In principle we can circumvent this parallax problem and construct maps with resolution 
comparable to the beam width of the experiment by assuming that the motion of the TNO on 
the Celestial sphere is completely due to the motion of the experiment and then accounting 
for the shift between the extragalactic CMB and the TNOs. This is extremely complicated 
but possible for low instrument noise and enough interlocking observations. The details will 
be presented in a future paper. 

When we change the distance to a group of TNOs from D to D + x, the fractional 
change in the TNO contribution to the observed CMB intensity is 

51 _ D 

T ~ 2x + xyD' ^^^^ 

The first-order term is a good approximation for D ^ x, but for shorter distances the second- 
order term will also contribute. For D ^ x we will be able to produce a two-dimensional 



^If the beam were completely filled then the distance modulation would not change the observed surface 
brightness. 
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map of the optical depth from small objects using this method, but for distances such that 
the second-order term is significant we will also be able to determine the correct distances 
to these objects and produce a three dimensional map. The number of TNOs in the beam, 
which is also unknown, is completely degenerate with the distance D when D ^ x and it is 
the second order variation that breaks this degeneracy and allows us to reconstruct both the 
beam filling fraction and the distance to the TNOs. We consider the potential of this effect 
to constrain TNO properties in the next subsection. 



4.1.4- Instrument Noise 

Each pixel will have a fluctuation due to instrument noise since the system temperature 
is above the CMB temperature. For reference we will describe instru ment noise in tempera- 
ture u nits (ur) and relate it to uncertainty in the observed intensity as (ITegmark fc Efstathiou 



19961) 



aB{iy) = o-T. (15) 



4.2. Potential Constraints 

We will now determine how well this technique can produce a map of TNOs for two 
model experiments, one with WMAP instrumental parameters and one with next genera- 
tion parameters. We will ignore the Doppler and Sachs- Wolfe effects discussed in §4.1.11 
and assume that their amplitudes can be determined with sufficient precision that we can 
correct for them. The Fisher matrix formalism will be employed to estimate these potential 
constraints (see Tegmark et al. 1997 for an overview of the Fisher Matrix formalism). 

The intensity in a given direction depends on our orbital position because of the portion 
coming from the outer Solar System. We can isolate this distance-dependent piece by dif- 
ferencing the observed CMB intensities in a given direction as measured at different orbital 
positions: 

AI^{n) = ^I^{h,Xi) - I^{h,Xj), (16) 
= ^[r(n,fi) - r(n,fj)][fi,,(TTNo) - 5i,(fcmb)], 
where the optical depth in a given direction is modeled as having an unknown constant 
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amplitude tq in addition to the unknown variable distance 

Note that Ttno(^) ^'Iso depends on D{n) through the thermal balance of Solar heating, Eq. 



In order to estimate the precision with which we can measure the distance to an trans- 
Neptunian object we will perform a Fisher matrix analysis. The Fisher matrix is defined 

where Pa = {To(n), D{n)}. Since Tq is unknown, we marginalize over this parameter, and find 
the constraints on D. Uncertainties in the CMB detector position and velocity, as well as 
uncertainties in the positions of planets within the Solar System, will all increase the overall 
uncertainty. We expect instrument noise to be the dominant source of uncertainty and will 
assume the noise in Eq. (|T8l) solely comes from the instrument noise. 



These limits allow us to detect the presence of TNOs in a given direction. To first order 
in Ax/D ~ 1AU/-D we are only sensitive to tq/D^, so our results would be completely 
degenerate between these two unknown parameters. This degeneracy is broken at second 
order in Ax/D. We can either try to make a two-dimensional map and constrain the TNO 
optical depth in each pixel or we can use second order variations in Ax/D to make a three 
dimensional map and determine both the optical depth in a pixel as well as the mean distance 
to TNOs in that pixel. Table [1] shows the potential constraints on tq and D that can be 
produced by a CMB experiment with a given pixel noise and that observes in ten frequency 
bands and measures a given pixel on the sky with ten different projected distances between 
-0.5 AU < Ax < 0.5 AU. 

The upper portion of the table corresponds to the Kuiper Belt. The plausibility of 
this technique strongly depends on the value of a. Our results imply that a map of the 
Kuiper Belt can be made if a ~ 2. In fact, there is a transitional value of amass > 5/3 
{ciradius > 1.22) corresponding to the point where the integral in Eq. ([2]) becomes dominated 
by the low mass end of the distribution. The lower portion of the table corresponds to the 
inner Oort Cloud. Due to the strong distance dependence of the effect producing a map of 
the inner Oort Cloud will be much more challenging. Assuming a = 1.9 and ctt = l/xK, 
there needs to be Mtotai = 0.185 x Msarth in order to produce a two dimensional map with 
signal-to-noise unity. 



- 19 - 



5. Discussion 

In this paper we describe two new tecliniques to constrain tlie cumulative mass contained 
in the outer Solar System. The distribution of mass in this portion of the Solar System may 
provide clues to both the early nature of and the formation scenario for the Solar System. 
Current techniques for finding these objects that rely upon the detection of reflected sunlight 
are only sensitive to the largest and closest members of this population, objects that are rare 
and therefore are likely to be special cases. Our new techniques directly probe the low mass 
end of the distribution that is otherwise quite difficult to detect. 

It has been suggested that observations of occultation events may be used to constrain 
the low mass end of the TNO distribution. Bodies in the Kuiper Belt can be detected when 
they transit in front of distant luminous objects, resulting in a brief dip in the observed 



brightness. The Taiwan- America Occultation Survey (TAOS) (iLehner et al.ll2006l ) seeks to 
conduct a census of large Kuiper Belt objects using background stars. The cadence of ob- 
servation required to detect the smallest objects presents significant technological challenges 
because of optical CCD read-noise. This same occultation m ethod has been ap plied to Scor- 



pius X-1, the brightest X-ray source near the ecliptic plane (jChang et al.ll2006l ). since X-ray 
detectors can be operated at a higher cadence than optical CCDs. However, the small num- 
ber of bright X-ray sources near the ecliptic limits the total volume of the Kuiper Belt that 
can be surveyed. 

Our method does not suffer from these limitations as the CMB detector technology 
is proven and the CMB can be accurately observed over nearly the entire sky (except for 
the Galactic plane). The main drawback of our proposed techniques is their reliance on 
the unknown TNO mass function. We must also be careful when interpreting our results 
because we are basically constraining the low mass end of the TNO distribution and then 
using those results to infer the total mass in the outer Solar System. If the slope of the 
high mass end of the mass function is steep (/? > 2), then most of the mass lies in the low 
mass end and this inference is valid. However, if the slope is relatively shallow (/? < 2) then 
we should only interpret our results as pertaini ng to the mass in sm all objects. Current 



observational results allow for both possibilities (iBernstein et al.ll2004l ) and hopefully future 
observations will determine the true slope. These concerns are in addition to those expressed 
in §3] regarding the possibility that rare high mass objects could cause the true total mass to 
differ from that inferred by extrapolating the mass function. Fortunately, the large objects 
are the easiest to optically detect and our proposed techniques can work together with these 
surveys to provide us with an accurate description of the outer Solar System. 

The next stage is to use the data provided by WMAP and analyze it as discussed in 
§1] in order to produce a map of Kuiper Belt objects. The data analysis, which must be 
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done with the time-ordered data, is a lengthy process, which is why it was not included 
in this paper. The potential to better understand the formation of our Solar System, and 
subsequently other planetary systems as well, will certainly justify the effort. 
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Mass 


Distance 


a 








IMe 


40 AU 


1.5 


\liK 


19300 


6470 


IMe 


40 AU 


1.66 


IfiK 


3280 


1100 


IMe 


40 AU 


1.8 


IfxK 


34.6 


11.6 


IMe 


40 AU 


1.9 


l^K 


0.575 


0.193 




10 AU 


2 


1//A' 


0.0192 


0.00647 


IMe 


40 AU 


1.5 




578000 


194000 


lAfe 


40 AU 


1.66 




98500 


33100 


IMe 


40 AU 


1.8 


SOfiK 


1040 


349 


IMe 


40 AU 


1.9 




17.2 


5.80 


IMe 


40 AU 


2 




0.577 


0.194 


50Me 


1000 AU 


2 


IfiK 


30.0 


10.0 


50Me 


5000 AU 


2 


lliK 


8700 


2910 



Table 1: Detection limits of the optical depth, ro(n), and distance, D(n), for Kuiper Belt 
and Oort Cloud objects. The optical depth and distance are partially degenerate, so the 
uncertainties presented here have marginalize over the second degenerate parameter. In 
order to make a meaningful map of the outer Solar System the fractional uncertainties must 
be less than one. 



